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We derive a Hamiltonian for a two-leg ladder which includes an arbitrary number of charge and 
spin interactions. To illustrate this Hamiltonian we consider two examples and use a renormalization 
group technique to evaluate the ground state phases. The first example is a two-leg ladder with 
zigzagged legs. We find that increasing the number of interactions in such a two-leg ladder may 
result in a richer phase diagram, particularly at half-filling where a few exotic phases are possible 
when the number of interactions are large and the angle of the zigzag is small. In the second 
example we determine under which conditions a two-leg ladder at quarter-filling is able to support 
a Tomanaga-Luttinger liquid phase. We show that this is only possible when the spin interactions 
across the rungs are ferromagnetic. In both examples we focus on lithium purple bronze, a two-leg 
ladder with zigzagged legs which is though to support a Tomanaga-Luttinger liquid phase. 
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I. INTRODUCTION 

Ladder systems are well known for their many novel 
properties and their relative simplicity makes them an 
ideal candidate for much theoretical workfii^ Several ex- 
perimental systems are known to be of or dominated by a 
ladder-type structure, and theoretical studies have been 
able to make reasonable predictions about the phases, 
symmetries and transport properties of these materi- 
alsi^i^i^ A common procedure used to solve ladder sys- 
tems is a perturbative renormalization group (RG) treat- 
ment, followed by a non-perturbative bosonization of the 
relevant interactions. The combination of these two com- 
plimentary techniques allows one to go beyond the usual 
mean-field approaches when determining ground states 
and excitations in the low-energy regime.^ Some studies 
using these techniques have revealed exotic phases, such 
as a staggered-flux phase'' and a resonant-valence-bond 
liquidi^ 

In a recent experiment,— it was demonstrated that 
Tomanaga-Luttinger liquid (TLL) behavior appears in 
lithium purple bronze Lio.gMoeOir (LPB). It is rather re- 
markable that typical TLL scaling appears to exist over 
a wide range of temperatures. In this reference it was 
claimed that RG flows quantitatively reproduce the ex- 
perimental data, but the bare interactions which lead to 
this solution were not discussed. This exciting develop- 
ment inspired us to revisit the well-known two-leg ladder 
system, modified to describe a realistic interaction pro- 
file while also taking into account the ladder geometry, 
as shown in Fig. [T] 

A standard two- leg ladder is shown in Fig. [T]Ja). The 
hopping strengths between nearest neighbors on the same 
leg and nearest neighbors on the same rung are t and t±_ 
respectively. For on-site interactions the charge and spin 
interactions take the same form and can be described by 
a single parameter U. In general, interactions between 
two different lattice sites along the same leg are described 
by X||„, while interactions between two lattice sites on 
opposite legs are described by X_l„. The charge and spin 




FIG. 1: (a) A standard two-leg ladder with hopping strengths 
t and t± and several electron-electron interactions defined 
by U and X where X — V, J . (b) A zigzag two-leg ladder 
where the legs are bent to make a constant angle (f). Electron- 
electron interactions can be defined similarly to the standard 
case. 



interactions are represented hy X = V, J respectively 
and the integer index n describes the rung difference be- 
tween the two sites. Therefore, any set of generic short- 
range or quasi-long-range interactions can be described 
by the bare interactions U, X|j„ and X±n- Although 
there have been extensive theoretical investigations on 
electronic correlations in two-leg ladders f^iii^iiSiiiiii^ most 
of these studies only consider nearest-neighbor (or next- 
nearest- neighbor) interactions. 

The standard two-leg ladder lies in a two-dimensional 
plane, but there are a number of experimental systems 
which contain a two-leg ladder which is warped in some 
fashion. ^^'^"^ We consider a ladder which has been com- 
pressed so that the legs form a zigzag with a constant 
angle 0, as shown in Fig. [TJb). One example of such 
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a lattice is LPB which has 4> ^ 7r/2ii^ By including the 
geometric structure of the two-leg ladder we have an ad- 
ditional variable (j> with which to investigate the ground- 
state phase diagram. It is easy to see that for extremely 
short-range interactions the zigzag angle does not play 
any significant role since the on-site interaction U domi- 
nates. However, in ladder materials the interaction is of- 
ten expected to be quasi-long-ranged and in these cases 
is important. With the inclusion of this zigzag lattice ge- 
ometry, as well as the generic interactions we hope to not 
only be able to fully describe a TLL phase but to also 
discover other exotic phases such as an /-density wave 
and a staggered-flux phase. 

The Tomanaga-Luttinger liquid is a special case 
amongst all the possible phases of a two-leg ladder. 
In sharp contrast to ordinary Fermi liquids, electron- 
electron interactions in TLL cause the single-particle ex- 
citations (the so-called quasi-particles) to become unsta- 
ble. Instead one finds bosonic spin and charge excitations 
which propagate independently of each other and with 
different velocities, a phenomenon known as spin-charge 
separation. Many theoretical and experimental studies 
have discussed TLL phase in several different ID or quasi- 
ID systems such as weakly-coupled chains or wires,— 
carbon nanotubeSj and the edges of two-dimensional 
systemsii^iSS With screened charge interactions theoret- 
ical studies have shown that TLL are generally expected 
in odd-leg ladders^i but not in even-leg ladders except 
under unphysical conditions such as attractive interac- 
tionsJi These general trends make the TLL scaling be- 
havior observed in LPB^^ a little unexpected, though 
as the two legs in LPB are almost independent {t ^ t±) 
it is certainly not impossible. There are two plausible sce- 
narios for the observed TLL- like behavior in LPB. The 
first scenario is that the ground state is a true TLL and 
that the interaction profile and the ladder geometry in 
LPB result in an unusual set of bare couplings that flow 
towards the TLL phase under RG transformations. The 
alternative possibility is that the ground state is not a 
TLL but closely resembles a TLL over a wide range of 
temperatures. In an attempt to solve this puzzle we will 
use LPB as an example when determining the phases of 
our two-leg ladder with zigzagged legs and generic inter- 
actions. 

This paper is organized in the following way: In Sec. 
II, we introduce the two-leg model that contains various 
charge and spin interactions. Starting from the lattice 
model, we briefly describe the chiral decomposition, cur- 
rent algebra, computation of initial couplings and the 
bosonization. In Sec. Ill, we generalize the theoretical 
approach to the two-leg ladder with general zigzag angles. 
We also introduce different order parameters to charac- 
terize the ground states. The complimentary combina- 
tion of the RG method and the bosonization techniques 
allows us to obtain the phase diagrams for different in- 
teraction profiles and bending angles. We consider two 
examples, the half-filled case with t = t± and the quarter- 
filled case with t t±, with the latter case corresponding 



to LPB. In Sec. IV, we make use of the general theoret- 
ical framework developed in previous sections and try to 
determine an appropriate interaction profile for a TLL in 
LPB. We perform detailed and extensive numerical anal- 
ysis and compute the temperature-dependent TLL ex- 
ponent. Finally, we conclude our numerical results and 
discuss their connections to experiments. 

II. MODEL 

We consider a two-leg ladder with quasi-long-range 
charge and spin interactions. The Hamiltonian contains 
non-interacting hopping as well as charge and spin inter- 
actions over different ranges. Thus, it is natural to divide 
the Hamiltonian into six parts, 

H = Ho + Hu + Hy^ + + Hj^ + Hj^^ . (1) 

The first term Hq describes hopping along the legs of the 
ladder with hopping strength t, and along the rungs with 
hopping strength t^, 

Ho = - t^ic],^Cj^i+j^)„ + h.c.) 

jla 

-ti-Y.^c\ia^iia + h.c.). (2) 

In 

The subscript of the fermion operator Cjic describes leg 
number j = 1, 2, rung number /, and spin a =1, |. 

For the on-site interaction, the difference between the 
charge and the spin parts vanishes and 

Hu = uJ2^jimn (3) 

jl 

where rijia- — c^ji„Cjia and U is the interaction strength. 
Now we classify the more general charge interactions. 
Many theoretical studies consider perpendicular nearest 
neighbor interactions across single rungs, i.e., between 
sites (j, /) and (j, /) where j denotes the opposite leg of j, 
as well as parallel nearest neighbor interactions between 
neighboring sites on the same leg, i.e., between sites (j, I) 
and (j, Z ± 1). A few studies also consider next-nearest 
neighbor interactions which act diagonally across one pla- 
quette, i.e, between sites and (j, Z ± 1). Here we 

consider all charge interactions between sites (j, Z) and 
(j',/') for which \l - V\ < N. Note that we have in- 
troduced a "hard" cutoff length for the interaction 
profile. The perpendicular Hamiltonian describes inter- 
actions between sites on different legs 

N 

Hv^ ^ X! X! '^^nnjlanj(i+ri)a' (4) 

where Vj^„ is the interaction strength between sites (j, I) 
and (j, I + n). The parallel Hamiltonian describes inter- 
actions between sites on the same leg 

N 

Hv^\ = X] 51 ^nnjlanj(l+n)a' (5) 
n=l jlaa' 
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where Vj|„ is the interaction strength between sites (j, /) 
and (j, I + n). 

Following the same classification the spin interactions 
are contained in two parts, Hj^^ and Hj^^ . Like iJy^ and 
i?yii , the spin interaction Hamiltonians describe interac- 
tions between any two sites which are N or less rung 
positions distant from each other. For spin interactions 
between sites on different rungs, 

N 

n=0 jl 

where Jl„ is the interaction strength between sites (j, /) 
and (j, I + n) and the spin operators are 

aa' 

where t — {Tx,Ty,Tz) are Pauli matrices. For spin inter- 
action between different sites on the same leg, 

JV 

-H",7|| = XI X! -^IW^Ol- ■ ^j(l+n) (8) 
n=l jlaa' 

where J||„ is the interaction strength between sites (j, /) 
and (j, I + n). 

We now follow a standard procedure which involves 
decomposing the lattice fermions into pairs of chiral 
fermions with linear dispersion. As this procedure is 
well explained elsewhere^ we will only give a brief ex- 
planation. Firstly, the hopping part of the Hamiltonian 
Hq is diagonalized into a bonding and antibonding band, 
dqia = [c2ia - (-l)*ci;o^]/V2 with q = 1,2, then after a 
Fourier transform we can determine the band structure 
Eg = (— l)'tj^ — 2tco&kq as a function of momentum kq. 
The Fermi momentum kpq = cos~^[{—iJ,+ {—\)H±)/2t] is 
uniquely determined by the chemical potential /u. As we 
are only interested in the low-energy behavior the fermion 
operators which diagonalize the hopping Hamiltonian can 
be linearized about the Fermi point by introducing chiral 
fermion fields, aqia ^ ipRqiae''^''''^ + V'i.j/cre"*''^''. Taking 
the continuous limit of the discrete lattice index I, we can 
define the Hamiltonian density Ti from H ~ J dVH. The 
hopping part of the Hamiltonian density in terms of the 
chiral fields is rather simple, 

Ho = - XI ^li^Rqladl^Rql'y - ^IqlA'^Lqla) (9) 

go- 
where the Fermi velocity is Vq — dEq/dkq at kq — kpq. 

The interaction part of the Hamiltonian density 7i/ = 
Tiu + 'Hv± + Ti-Vn ~^ ^■J± + ^J\i be expressed in terms 
of the currents 

Jpqq' ^Wpqa^Pq'a, 3 Pqq' = ^i^lqa'^ aa' i^Pq' a' 

Ipqq' —^^JPqa^aa'i^Pq'a' , ^Pqq' — jtp Pqa{(-f) aa' Pq' a> 

(10) 



where P = R, L. The antisymmetric matrix e is defined 
by = —£21 = 1 and £11 = £22 = 0. Each term in Hi 
is a product of two currents so that the Hamiltonian is a 
function of four-fermion interactions, 

H/ =bqqiJRqq'JLqq' " b'^q' ^ Rqq' Lqq' 

~^ fqq"^Rqq'^Lq'q' ^ fqq'-^Rqq-'^Lq'q' 

+ Uqq'llq'lLqq' " W-lLq,'- (H) 

The couplings of the four-fermion interactions bqqi , fqqi 
and Uqq' define the scattering amplitudes between bands 
q and q' . Backward scattering is represented by b and 
from a gradient expansion the bare coupling strength can 
be shown to be 

N 

b'qq =U + y±o + 2 XlC^I" + ■^^")(2 " cos 2nkFq)] 

n=l 
AT 

- |^_L0 - I ^[i.J\\n + J^n)COS2nkFq] 

n=l 
N 

bqq =f/ + Fio + 2 X ii^n + V±n) COS 2nkFq] 
n=l 
N 

- jJ±0 -J2^iJ\\n + J_Ln)(l + ^C0s2nkFq)] 

Tl=l 

N 

bi2 =U - V±Q + 2^[{V\in - yj.„)(2cosnfc_ - cosnfc+)] 

n=l 

JV 

+ I-/-LO - I ^[{J\\n - J±n)cOSnk+] 
n = l 
N 

bi2 =U - Vi_a + 2 XK^I" - ^^n) cosnfc+] 

n=l 

N 

+ I'/io - ^[{J\\n - J_L„)(cosnfc_ + icos?ifc+)] 

(12) 

where k± = kpi ± fc_F2- The symmetry of the system 
requires &12 = ^21 and at half-filling fi — so kpi + fcF2 = 
TT which sets bn — &22- Forward scattering is represented 
by / with the bare coupling strength, 

N 

fP^ ^[/ + 2 ^[V|l„(2 - cos nk+) + Vl„(2 -I- cos nk+)] 

n=l 

N 

+ 'iV±o + |Jio - ^^[{J\\n - J±n) COSnfc+] 
n=l 

N 

=U - V^o + 2Y,[iV\\n - V^n) cosnk+] - i J^o 

n=l 

AT 

— Xt'^ll"^"'" ^ 5 cosnfc+) — J±n(l — \ cosnfc+)]. 

(13) 
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Symmetry requires /12 = /21 and in order to avoid dou- 
ble counting we set fqq — 0. Umklapp scattering is 
represented by u and only present at half-filling where 

kpi + kp2 ~ TT, 

N 

<i =C/-yio + 2^[(l^|„-yi„)(-l)"] 

n=l 
N 

n=l 

N 

<2 =2C/ + 2^[V||„(cos2nfci.i + (-f)") 

n=l 

+ V^n{coB2nkFi - {-IT)] 

N 
n=l 

- Ji„(cos2nfcFi - (-1)")] 

N 

<2 =2yio + 2^[l^|„(cos2nfcj.i- (-1)") 

n=l 

+ V±nicoa2nkFi + (-1)") + ^J±o 

TV 

+ kJ2^J\\n{cos2nkFi-i-ir) 

n=l 

+ Ji„(cos2nfcFi + (-1)")] (14) 

with U12 — U21 and un ~ U22 from symmetry. At half- 
filling, the particle-hole symmetry ensures we have nine 
unique coupling constants, b^^ = 6221 ^11 = ^22 j ^121 ^12 ' 
fi2, fi2i '^11 = "221 "12 ^'^'i ""12- Away from half- filling 
the Umklapp interactions vanish but we no longer have 
611 = ^22 so we have eight different coupling constants. 

The Hamiltonian is more easily analyzed if the chiral 
fermion operators are replaced with boson operatorsi^^ 
The bosonized fields and <y9^± with v = p,a represent 
a variety of quantum numbers. The subscript represents 
total {+) or relative (-) charges or spins (p or a respec- 
tively) between the two bands while is a displacement 
field and is a phase field. The total bosonized Hamil- 
tonian density H = Ho + Hi is 



on 

~ 26^2 cos (pp_ cos 9cr+ + 2 cos 9cr+ {b^i cos 9cr- 
+ fi2 COS (/3o-_ ) — cos (y9p_ cos ^cr- + ^12 COS fcr-) 

— 2u1^ COS 6p^ COS (/5p_ — 2m^^2 cos 6p-^ cos 9^^ 

— COS 0p+('U^2 cosher- + M]^2 COS ((9ct-) (15) 

where = ^12 ^ ^12 ^^"^ '^12 ~ ""12 ^ "12 ■ '^^^ Luttinger 
parameters and the Fermi velocities for the total/relative 



charge and spin sectors are 

^ ^ l 27T{v,+V2)-m,+b^22)/2±f^2\ .... 

^ 27r{v,+V2) + m,+b'^,)/2±f-,] ^ > 

v.± =^47r2(«i + V2r - [ib'(, + W22)I2 ± rx2?l^n. 

Note that the highly symmetric bosonized form in Eq. 
([T5)) is possible only for degenerate velocity v\ = V2- This 
is always true at half- filling but not at generic fillings. 
The LPB two-leg ladder we are interested in is at quarter- 
filling,-^ but as it consists of nearly independent chains 
with vanishingly small inter-chain hopping t±_ <ti t, the 
Fermi velocities are nearly degenerate vi ^ V2 and the 
above bosonized Hamiltonian is valid. 

The RG flow equations of the couplings are of the 
form dgi/d£ — X^jfc M'^ djdt^^ where A!'^ is a constant ten- 
sor that can be computed from operator product expan- 
sions. ^'^■^ All RG flow equations are solved simultaneously 
with the initial conditions at ^ = given in Eqs. (|12llf 4p . 
The ground state phase is determined from a non-zero or- 
der parameter, such as electron density or current flow, 
which can be evaluated using the solutions of the RG flow 
equations. On substituting the RG solutions into the 
bosonized Hamiltonian Eq. (|f 5p the Hamiltonian may 
be minimized by a speciflc set of pinned bosonized flelds 
(while other fields remain free to adopt any value), thus 
defining the ground state. For example, if 6^2 flows to a 
non-zero value then 6^2 cos (^p_ cos 0^+ in Eq- (|15p could 
minimize the Hamiltonian by pinning (^p_ , 9„^ = mn for 
some integer m. In order to maintain this minimum, m 
can change by integral values, which describes an excita- 
tion over some finite energy gap. If instead 6J2 flows to 
zero the term bi2 cos ipp- cos 9^+ is irrelevant. Any field 
that remains unpinned when the Hamiltonian is mini- 
mized may describe a gapless excitation. The pinned 
fields may be substituted into the order parameter equa- 
tions, once they are appropriately bosonized, to deter- 
mine the phase. Therefore, each phase can essentially be 
defined by a unique set of pinned boson fields, which are 
directly related to the solutions of the RG flow equations. 
Note that it is only the coefficients of the sinusoidal terms 
which ultimately determine the gapped excitations and 
therefore and fi2 are the only couplings which can be 
non-zero in a fully gapless phase, i.e., a TLL. 



III. ZIGZAG TWO-LEG LADDER 

Now we would like to incorporate the realistic ladder 
geometry into the above theoretical model. We consider 
a two-leg ladder in which the legs zigzag parallel to each 
other, as shown in Fig. [Hb). The bent legs make a con- 
stant angle 4>. For unscreened charge interactions the 
interaction strengths between two sites are inversely pro- 
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portional to the distance between them so, 
X 



an sin 6/2 



n = 2,4,6,... 



X\\n — 



X 



X \ n — - 



a^l + (n2 - l)sin(/)/2' 
X 



X 



±n 



ay 6^ + v? sir? cf)/2 
X 

a^l + (52 + (n2 „ 1) sin 0/2 ' 



n = 1, 3, 5, . 
n = 0,2,4,... 



n = 1,3,5,, 



(17) 



where X — V, J , a \s the distance between neighboring 
lattice sites on the same leg and the distance between 
lattice sites on the same rung is a5. We shall assume the 
ladder consists of square plaquettes with 5 = 1. Gen- 
erally we would like interactions beyond the cutoff, i.e., 
with n > A^, to be less strong than interactions within 
the cutoff, but when (p is very small this may not be the 
case. This issue may be avoided by defining different cut- 
offs for interactions along a leg and interactions between 
legs, but as the small values of for which this problem 
occurs are quite likely not experimentally attainable we 
will continue to use just one cutoff N . 



A. Order parameters 

At half-filling the phase of a two-leg ladder could be 
one of four density wave phases or one of four Mott insu- 
lator phases i^iii We first discuss the density wave phases, 
the charge density wave (CDW), the staggered-flux (SF) 
phase, the p-density wave (PDW) and the /-density wave 
(FDW), and their associated order parameters. A CDW 
has a non-zero variation in the average electron density 
per site which is defined by 



riji 



(18) 



At half-filling the average electron density of a two-leg 
ladder is one electron per site, but in a CDW the sites are 
alternatively unoccupied or fully occupied by two elec- 
trons. To define current fiow we use 

CJ 
CJ 

a 

which describe currents along rungs between sites (j, V) 
and (j, Oj along legs between sites (j,^) and (j, ^ -I- 1) 
and along the diagonals of the plaquettes between sites 
(j, Z) and (j, Z -I- 1), respectively. If the first two currents 
are non-zero we have a SF phase which is characterized 



by alternative clockwise and anticlockwise current flows 
around plaquettes. If the third current is non-zero we 
have a FDW which is characterized by currents zigzag- 
ging across the diagonals of the plaquettes. Two kinetic 
order parameters are 

(7 



where B\\ji describes interactions along legs and Bdji 
describes interactions across diagonals. A PDW is de- 
fined by non-zero which implies dimerization be- 
tween neighboring sites on the same leg. While non-zero 
Bdji does not formally define any phase, it tends to be 
non-zero in a CDW and describes dimerization between 
sites of equal electron density. The kinetic energy across 
rungs 



(21) 



is always zero. 

Away from half-filling the situation is a little different 
with there being only two possible density waves phases. 
In this case a CDW (SF) and a PDW (FDW) coexist in 
a single phase, and for simplicity we name this phase a 
CDW (SF) phase. At half-fiUing the CDW (SF) and the 
PDW (FDW) only differ by the pinned value of the total 
charge displacement field 0p+ . Away from half-filling the 
Umklapp terms are removed and this provides an addi- 
tional symmetry, resulting in an unpinned 9p+. When 
dp^ is unpinned the CDW and PDW may coexist with 
the relevant order parameters, Uji and being simul- 
taneously non-zero. Similarly, the SF and the FDW may 
coexist and all three currents j_Lji, and jdji will be 
simultaneously non-zero. 

If all the order parameters discussed above vanish then 
we may have a Mott insulator or a superconductor state. 
An s-wave superconducting order parameter can be de- 
fined by 



^slj - Cjl^Cjll ^ ^Pql 
^ P<1 



(22) 



where Ap,; = ipPg^i^Pqii ^® pairing operator of the 
chiral fields. The d-wave order parameter across the 
rungs is 



Ad±i = cu^C2ii - i^(-l)«+iAp,,. 



(23) 



Pq 



As the names imply, Agij is non-zero in an s-wave super- 
conductor (S-SC) while Ad±i is non-zero in an d-wave su- 
perconductor (D-SC). On bosonizing the superconduct- 
ing order parameters it can be seen that they can only be 
non-zero away from half-filling where the boson field 0p-|_ 
is unpinned. At half-filling the total charge displacement 
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is pinned and both A^ij and Ad±i vanish, and provided 
all previously discussed order parameters are also zero 
we may have a Mott insulator. The Mott insulator at 
half-filling is defined by non-zero Apqi and, like a super- 
conductor, is defined in terms of a pairing symmetry. If 
Afli;Aj^2! > we define the Mott insulator as s-wave, 
but if Afli/Ajj2i < we define it as d-wave. Two types 
of s-wave and d-wave Mott insulators exist, one with 0p_|_ 
pinned to an even multiple of tt, named S-Mott and D- 
Mott, and the other with 9p+ pinned to an odd multiple, 
named S'-Mott and D'-Mott. A difference in 9p+ of tt 
represents a half-plaquette shift in the center of mass of 
the paired chiral fields with the D and S-Mott pairing 
being across rungs and the D' and S'-Mott pairing being 
across the diagonals of the plaquettes. 

There are a few other possible phases in the two-leg 
ladder. For example, the phases transitions between any 
two phases may be though of as phases in their own right, 
but as they only exist over a vanishingly small parameter 
range we will not discuss them here. Finally we should 
mention the TLL, although we shall discuss this phase 
in more detail in Sec. IIVI In a TLL all bosonic fields 
arc unpinned and because of this all order parameters 
discussed above are undefined. 



B. Half- filling 



Our first example of the zigzag two-leg ladder is the 
case of half-filling fi ~ with equal leg and rung hopping 
t = t±. Phase diagrams constructed from the solutions of 
the RG flow equations are given in Fig. [2] and clearly both 
(j) and play a significant role. We always assume on-site 
interaction U = 2, although in the results presented here 
it is the ratios V/U and J/U which are important in de- 
termining the phase, rather than the actual values of U, 

V and J. When there are charge interactions V but 
no spin interactions J = increasing N from 2 to 10 will 
allow CDW and S-Mott phases to emerge while signifi- 
cantly reducing the range of the D' and S'-Mott states. 
For spin interactions J ^ but no charge interactions 

V ^ 0, when iV = 2 only a D-Mott and a PDW phase 
are possible, and only at quite small angles, < 7r/3. As 
A^ is increased to 10 one still requires (f) < Tr/3 to obtain 
anything but a D-Mott phase, but at these small angles 
several phases are possible. Of particular interest is the 
emergence of a FDW as this phase has not previously 
been predicted in a two-leg ladder at half-filling under 
any physically possible scenarios. Although the angle re- 
quired to obtain a FDW is quite small it may be possible 
to construct an appropriate lattice using cold atoms i?'* 
The phase diagram for A = 10 and = 7r/4 with vari- 
able V and J is shown in Fig. [3] Although the PDW 
dominates, there is quite a substantial FDW region. 
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FIG. 2: Phase diagrams at half- filling with t = = = 2 
and (a) iV = 2, J = (b) AT = 2, 1/ = (c) iV = 10, J = 
(d) = 10, = 0. Dashed lines indicate where is small 
enough so that some interactions not considered are larger 
than some which are considered. 
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17 = 2 and = 7r/4. 



C. Quarter- filling 

For our second example of the zigzag two- leg ladder we 
assume quarter- filling which sets \x = —^J2t'^ — t\ and 
we also assume t ^ 1 = 0.01. This case is designed 

to correspond to LPB when (f) ^ 7r/2. We again use U = 
2, although, as before, it is the ratios V/U and J/U which 
ultimately determine the phase. In Fig. |4] we present a 
number of phase diagrams. Despite the small SF phase, 
the J 7^ 0, I^ = case is not particularly interesting 
as extremely small values of are required if any phase 
but a D-SC is to be observed, particularly when N = 2. 
Unlike the half-filled case, increasing A^ from 2 to 10 does 
not cause new phases to emerge, although the SF phase 
does appear at a larger value of 4>- In the V ^ 0, J = 
case increasing A^ from 2 to 10 decreases the complexity 
of the phase diagram, causing the CDW phase to expand 
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FIG. 4: (Color online) Phase diagrams at quarter-filling with 
t = 1, t± ^ 0.01, (7 = 2 and (a) iV = 2, J = (b) iV = 2, 
V = (c) N = W, J ^ (d) N = W, V = 0. Dashed lines 
indicate where </> is small enough so that some interactions 
not considered are larger than some which are considered. 



and the S-SC phase to vanish. 

The region in Fig. [UJa) which is marked in red de- 
scribes a region of unusual scaling. The phase in this 
region is either D-SC or S-SC, with the D-SC to S-SC 
phase transition running approximately through the cen- 
ter. Each phase is characterized by a unique set of RG 
solutions of the eight coupling constants and generally, 
while renormalizing, the coupling constants flow gradu- 
ally towards this final final solution. In the red region the 
coupling constants do not initially flow towards either a 
D-SC or S-SC solution but instead towards a solution 
typical of the D-SC to S-SC phase transition. This scal- 
ing behavior continues as £ increases, but at some point 
there is a sudden change and the RG will flow rapidly to 
either a D-SC or S-SC solution. This scaling behavior is 
typical when extremely close to a phase transition, but 
it is not usually observed in regions as large as the red 
region in Fig. Ufa). It is quite possible that this region 
could be mistaken for a TLL phase, as we shall explain 
in more detail in the next section. Note that this re- 
gion is very close to (/) = 7r/2 so it may explain the TLL 
observations in LPBi^ 



IV. TOMANAGA-LUTTINGER LIQUID 

In the previous section we constructed various phase 
diagrams while assuming physically realistic conditions, 
yet did not observe a TLL. In this section we look more 
closely at what is required for the RG equations to flow 
towards a TLL solution. Wc simplify the problem a lit- 
tle by considering a two-leg ladder with the interaction 
cutoff iV = 1. Note that in this limit, the zigzag an- 
gle (j> does not affect the initial values of the couplings 
and thus can be ignored. This model was discussed in 
Ref. d to describe LPB and solved using RG flow equa- 



tions equivalent to the ones used here. Note that one 
of the key features for TLL is the critical exponent a of 
the single-particle density of states. Quasi-particle exci- 
tations are not found in a zero-temperature TLL so the 
single-particle density of states p(e) at energy e should be 
suppressed near the Fermi energy ep. This suppression 
is expected to follow a power law p(e) cx |e — e^j" for 
some positive constant a as the temperature approaches 
zeroi^ 

In Ref. H, it was argued, both experimentally and the- 
oretically, that the nature of the critical exponent a in- 
dicates that LPB has a TLL phase. A remarkable agree- 
ment was found between the experimental value of a and 
the theoretical value obtained from the RG solutions, but 
what electron-electron interactions would provided the 
required initial conditions of the RG equations were not 
stated. Here we will discuss the electron-electron inter- 
actions which may support a TLL in a LPB-like two-leg 
ladder. 

We can evaluate a numerically at different tempera- 
tures from the coupled RG equations. Note that the 
the temperature scales as T = Tqc"^ under RG trans- 
formations, where Tq is the initial (bare) temperature. 
Therefore, when calculating the couplings' flow with the 
logarithmic length scale £, we can compute the critical 
exponent a at different temperatures. It is known that 
the exponent takes the form, 

a = {Kp+ -f K;1 + + - 4)/8. (24) 

If it approaches a constant during the RG analysis we 
have a hint of TLL behaviour.^ For convenience we sep- 
arate this critical exponent into two parts, a = a+ + 
where 

a± = {Kp± + - 2)/8. (25) 

As has been discussed previously, the couplings in front 
of the sinusoidal terms in Eq. (jisp determine the en- 
ergy gaps and thus the nature of the phase. If none of 
these couplings become relevant under RG transforma- 
tion, the ground state is gapless and is characterize by 
the so-called Luttinger parameters Kp± and Ka-± in the 
charge and spin sectors. In this case only the first line of 
Eq. remains, which corresponds to the TLL Hamil- 
tonian.^^ If at least one of the coefficients of the sinusoidal 
terms does not flow to zero we have any one of the Mott, 
SC, or density wave states discussed above. Three cou- 
plings, b^q with q = 1,2 and /fj, are not coefficients of 
sinusoidal terms so need not vanish in a TLL. From the 
RG equations it can be seen that these three couplings 
will remain roughly constant when, and only when, all 
the other gap-inducing couplings are irrelevant. ^'^^ Fur- 
thermore, only these three couplings appear in Eq. ([24]) 
which defines a. This is in agreement with what we have 
already stated about a TLL, i.e. the RG solution of a 
must flow to a constant value. 

We can make some comments about the general behav- 
ior of a. From the RG flow equations we determine that 
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always decreases, but /fj always increases. In fact, 



K2 



remains constant in RG flows. There- 
fore Kp+ must be a constant and Kp- always increases. 
This implies a+ is a constant and the flow of the expo- 
nent a is essentially determined by a_. The minimum 
of a will be when a = q;+, which corresponds to a_ = 
and Kp^ = 1, or equivalcntly b^^ + = 2/^2. When 
Kp^ < 1 (or + 622 > S/fj) ct^ will decrease, as must 
a, but when Kp^ > 1 both a- and a will increase. The 
point b^i + 622 = 2/^2 is more significant that just the 
turning point of a. The RG flow equations indicate that 
once this point has been reached our system cannot be a 
TLL and a will increase at an increasing rate. So, once 
Ki + ^22 < 2/f2, or equivalcntly Kp^ > 1, in the RG 
flows, the TLL phase is unstable and some energy gaps 
will appear. However, it is important to emphasize that 
with b^^ + 622 > 2/C2 satisfied we may have a TLL but it 
is not guaranteed. 

Returning to our specific example of LPB, we again 
assume quarter-filling and set t = 1, t± — O.OI. In this 
case kpi, kp2 ^ n/A. Using the initial conditions in Eqs. 
(fT2| and (fT3|) it can be shown that the condition b^^ + 
622 > 2/C2 is equivalent to 



N/2 
n=l 



N/2 

-1)"J±2„ < -|[V±o + 2^(-1)"Kl2„]. 



In order to compare our results with Rcf. 
consider the charge interactions V± = V±0; 



(26) 
8 we only 

Vd = V^i, 



Vj| — Vj|i and the spin interactions J± = J±o, and 
J|| = J||i which corresponds to a cutoff iV — 1, and so 
the necessary but not sufficient condition for a TLL re- 
duces to Jj_ < — We wish to restrict ourselves to 
physically possible cases so we must have V± > and 
therefore the spin interaction across rungs J± must be 
negative, implying ferromagnetic exchange coupling. In 
Fig. Owe present four phase diagrams, all of which show 
a LL may be obtained for negative J±. In all cases the 
condition J± < —^V± is satisfied when we have a TLL, 
but clearly it does not imply that we must have a TLL. 

If we wish to choose initial conditions which will en- 
able the RG flow of a to closely resemble the experi- 
mental data in Ref. @ then even greater restrictions 
are placed on J±. As the interaction strengths are un- 
known we must attempt to make a reasonable guess. 
We set [/ = 2 and to make the numerical search prac- 
tical we set V±_ ~ V\\. Then, we fit the experimental 
data by varying the bare values of the charge interac- 
tions V±,Vd and the spin interactions J^, Jy. Accord- 
ing to experimental dataS 0.6 < a < 1, with a = 1 
corresponding to the highest temperature measurement. 
So, we set a+ = 0.6 as this marks the minimum of 
a and this sets Kp+ = 0.15. The maximum value is 
set to a{To) = 1 so a^{To) = a{To) - a+ = 0.4 and 
Kp-{£ = 0) = 0.20. These conditions determine the ini- 
tial values of b^^ + 622 ± ^fi2 which in turn determine 
y^iVd for a given choice of J±,J\\. The RG flow of a 
which corresponds to the experimental data is shown in 
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FIG. 5: Phase diagrams at quarter-filling with t = \, = 
0.01, i7 = 2 and (a) V\\=Vi = J\\ = 0, (b) = Vd = Jy = 0, 
(c) l^x = 1^1 = J|| = 0, (d) Vl = = Vd = 
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FIG. 6: (Color online) The critical exponent a at quarter- 
filling with t = \,t^= 0.01, V = 2,Vj| = Vl and Jy = (1-f 2j) 
with j = 0, 1, 2. While the arrow indicates increasing Jy for 
the curves show here, it does not indicate a general trend. 
The lowest red line and the lowest green line diverge so are 
not TLL, while the other four lines approach constant values 
and imply a TLL state. 



Fig. m with the initial temperate Tq = 300 K. We find 
that a TLL with 0.6 < a < 1 is only obtained when J_l 
is quite large (significantly larger than U) and negative, 
which is unrealistic for LPB. This does not mean that 
a TLL phase is impossible for LPB as we must bear in 
mind that these simple one-loop RG solutions are not 
expected to be quantitatively correct and should really 
only be used for qualitative analysis. Consequently, forc- 
ing the theoretical value of a to fit the experimental data 
is not recommended and will not give a good prediction 
of the interactions in the lattice. 

Despite the limitations of this RG method it has had 
some success in predicting phases of various systems. 
Rather than attempt to obtain the exact experimental 
values of a one could simply try to replicate the line 
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shape of a as the temperature decreases. In Fig. [Tlja) we 
show that the general Hne shape observed in experiments 
is obtainable when Jj_ is not particularly large, although 
it must be ferromagnetic because we are still bound by 
the condition Jj_ < — |V_l if we wish to have a TLL. In 
Fig. [TJ^b) we show that when we do not have a TLL a 
may still adopt a variety of line shapes, some of which 
strongly resemble a TLL as their turning point is very 
close to T = 0, in particular the J\\ = —J± = 1/2 case. 
Also, by rescaling the interaction strengths it is possible 
to rescale almost any a to have a very low temperature 
turning point. 

If we rescale all interactions by the same factor R so 

that (t/, J„||, J,i±,K||,K_l) ~* {U,Jn\\,Jn±,Vn\\,Vn±)/R 

then, because the initial couplings are linear in the inter- 
actions we can define a new set of couplings gt in terms 
of these rescaled interactions, gi{i = 0) = Rgi{i = 0). 
Rescaling the RG flow equations in the same way gives 
dcji/di = '^jk-^i'^Qjf/k where £ = Ri and because it is 
the ratios (J„||, Jnj_, Vn±)/U which essentially deter- 
mine the phase both gi and gi should flow towards the 
same solution and eventually describe the same phase. 
However, this does not mean they will have the same 
scaling. The rescaled temperature is T = Tqc"' and so 
T = To(T/To)~^ and therefore, by choosing an appro- 
priate R we may rescale a so that its turning point is 
very close to T = and the phase may closely resemble 
a TLL over a large temperature scale. For example, the 
J|| = — = 4 curve in Fig. [7^b) has a turning point at 
r=115 K, but if we choose R = 2 we rescale to U = 1 
and J|| = — Jj^ = 2 which rescales the turning point of a 
to T = 44 K. Similarly, if we choose i? = 4 we obtain a 
turning point at T = 6.4 K. 



In Fig ([5]) we show the RG flow of the couplings for the 
U = 2, J|| = —J± = 1/2 case. These couplings mostly 
behave very much like one would expect in a TLL, with 
f^2j ^12 ^11 approaching zero while fejj^ and are 
fairly constant, resulting in a fairly constant a over a 
large temperature range. Only does not have typi- 
cal TLL behavior as it does not approach zero. Because 
of this the RG eventually flows away from typical TLL 
behavior and the couplings become large, in this case 
flowing towards a typical D-SC solution. In the previous 
section it was mentioned that the region outlined in red 
in Fig. HJa) is not a TLL, but may be mistaken for one. 
This is because these couplings scale similarly to the ones 
shown in Fig ([5]), specifically b^i and /fj (and therefore 
a) remain fairly constant over a significant £ range but 
at some point they make a rapid change and approach 
values typical of a superconductor. As the experimental 
data of LPB also hints at an increase in the critical ex- 
ponent a near T = 0,— it cannot be ruled out that the 
observed scaling is indeed a close crossover from TLL- like 
behavior to some superconducting or density wave phase 
near zero temperature. 
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FIG. 7: (Color online) The critical exponent a at quarter- 
filling with t = l,t±^ 0.01, f/ = 2 and (a) Jy = Vl = V|| = 
Vd^O resulting in a TLL, (b) Jy = ~J^, VI = V|| = = 
not resulting in a TLL. 
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FIG. 8: (Color online) The renormalized coupling constants 
(rescaled by 4nVq) at quarter-filling with t — 1, t± = 0.01, 
17 = 2, J|| = -Jx = 0.5 and Vj|,x,d = 0. 



V. CONCLUSIONS 

We have derived a Hamiltonian for a two-leg ladder 
which allows consideration of generic short-range charge 
and spin interactions. One can choose the interactions 
to extend only to nearest neighbors, or one can choose to 
have interactions which extend across several lattice sites. 
When increasing the range of the interactions the number 
of variables inevitably increases. Rather than considering 
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each interaction strength as an independent variable and 
deahng with all the associated problems, we can simply 
assume that the interactions are inversely proportional 
to the distance between lattice sites, thus keeping the 
number of parameters to a minimum. Using this method 
we were able to solve the RG solutions for any number 
of interactions while only needing three variables, C/, V 
and J, to describe the electron-electron interactions. 

The Hamiltonian derived here is applicable to several 
different materials, and not just those materials like LPB 
which have an obvious ladder structure. Carbon nan- 
otubes, for example, have an hexagonal lattice structure 
which may be mapped onto a two-leg ladder, and results 
obtained from two-leg ladder RG equations have been 
applied to carbon nanotubes with nearest-neighbor in- 
teractions.^^ However, carbon nanotubes are known to 
support long-range interactions so the Hamiltonian pre- 
sented here, with slight modifications, would provided a 
more accurate picture of the phases of a carbon nanotube. 

Our RG analysis of LPB is somewhat limited because 
we have no experimental data which gives any clear in- 



dication of the charge and spin interaction strengths. It 
is important to note that this data should not be ob- 
tained indirectly by attempting to fit the experimental 
flow of the critical exponent a to numerical solutions of 
a obtained from the RG equations as these numerical 
solutions are not expected to be quantitatively accurate. 
Because of these limitations we are unable to make a 
definite statement concerning a TLL phase in LPB. The 
observed behavior may be a true TLL phase, or it may 
simply be a non-TLL phase which strongly resembles a 
TLL over a large temperature range. The power of these 
one-loop RG solutions is that they are relatively sim- 
ple and tend to provide a qualitative description of the 
phase of the system. More quantitative accuracy may 
possibly be achieved from second-loop or higher order 
corrections;^ 
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from the National Center for Theoretical Sciences in Tai- 
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